Critical threshold and dynamics of a general rumor model on complex networks 
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We derive mean-field equations that describe the dynamics of a general model of rumor spreading 
on complex networks, and use analytical and numerical solutions of these equations to examine the 
threshold behavior and dynamics of the model on random graphs, uncorrelated scale-free networks 
and scale-free networks with assortative degree correlations. We show that in both homogeneous 
networks and random graphs the model exhibits a critical threshold in the rumor spreading rate 
below which a rumor cannot propagate in the system. In the case of scale-free networks, on the other 
hand, this threshold becomes vanishingly small in the limit of infinite system size. We find that the 
initial rate at which a rumor spreads is much higher in scale-free networks than in random graphs, 
and that the rate at which the spreading proceeds on scale-free networks is further increased when 
assortative degree correlations are introduced. The impact of degree correlations on the final fraction 
of nodes that ever hears a rumor, however, depends on the interplay between network topology and 
the rumor spreading rate. Our results show that scale-free networks are prone to the spreading of 
rumors, just as they are to the spreading of infections. They are relevant to the spreading dynamics 
of chain emails, viral advertising and large-scale data dissemination algorithms on the Internet, the 
so-called gossip protocols. 
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I. INTRODUCTION 

Many real-life technological, social and biological sys- 
tems have complex network-like structures, where ver- 
tices represent entities comprising the systems and edges 
the presence of some form of interaction or relationship 
between them. Some important examples include the 
Internet |U, 121 > the World Wide Web |3( , and social inter- 
action networks 0, H, El • ^ n recent years a significant 
amount of research has been devoted to empirical and 
theoretical studies of such networks 0, ^| . 

A large body of this work has been focused on finding 
statistical properties, such as degree distributions, aver- 
age shortest paths, and clustering coefficients that char- 
acterize the structure of complex networks. Another ac- 
tive area of research is on creating realistic models of such 
networks which can help us better understand the un- 
derlying organization mechanisms behind their formation 
and evolution. A third area of research, which is still at 
its early stages, investigates dynamical processes on com- 
plex networks with the aim of understan ding the im pac t 
of network topology on the dynamics 0, Il2l HH H3.ll5l| . 

A prominent example of such process is the spreading 
of rumors in social networks jlo . fl7L fl3 | . Rumors bear an 
immediate resemblance with epidemics and their spread- 
ing plays an important role in a variety of human affairs, 
such as stock markets politics and warfare. Rumor- 
like mechanisms also form the basis for the phenomena 



of viral marketing, where companies exploit social net- 
works of their customers on the Internet in order to pro- 
mote their products via the so-called 'word-of-email' and 
'word-of-web' pol l2l|. Finally, rumor-mongering forms 
the basis for an important class of communication proto- 
cols, called gossip algorithms, which are used for large- 
scale information dissemination on the Internet, and in 
peer-to-peer file sharing applications [2^, . 

Unlike epidemic spreading which has been extensively 
studied, quantitative models and investigations of rumor 
dynamics have been rather limited. Furthermore, most 
existing models either do not take into account the topol- 
ogy of the underlying social interaction networks along 
which rumors spread, or use highly simplified models of 
the topology. While such simple models may adequately 
describe the spreading process in small-scale social net- 
works, via the word-of-mouth, they become highly inad- 
equate when applied to the spreading of rumors in large 
social interaction networks, in particular those which are 
formed on the Web. Such networks, which include email 
networks |2^|, and social networking sites 0, typically 
number in tens of thousands to millions of nodes and their 
topology can show highly complex connectivity patterns 

013. 

In two previous papers 0, |2j| we investigated the 
dynamics of a classical model of rumor spreading, the 
so-called Daley- Kendall model [l7|[. on such networks. 
In this paper we describe a more general model of ru- 
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mor dynamics on networks, which unifies our previ- 
ous model with the Susceptible-Infected-Removed (SIR) 
model of epidemics. We formulate the dynamics of 
this model within the framework of Interacting Markov 
Chains (IMC) 26], and use this framework to derive 
a set of equations that describe on a mean-field level 
the dynamics of the model on complex Markovian net- 
works. Using analytical and numerical solutions of these 
equations we examine both the steady-state and the 
time-dependent behavior of the model on Erdos-Reyni 
(ER) random graphs, as a prototypical example of ho- 
mogeneous networks, and strongly heterogenous scale- 
free (SF) networks, both in the absence and presence of 
assortative degree correlations. 

Our results point out to several important differences 
in the dynamics of rumor spreading on the above net- 
works. In particular, we find the presence of a critical 
value for the rumor spreading rate (a rumor threshold) 
below which rumors cannot spread in random graphs, 
and its absence in scale-free networks with unbounded 
fluctuations in node degree distribution. We find that 
the initial spreading rate of a rumor is much higher on 
scale-free networks as compared to random graphs, and 
that this spreading rate is further increased when assor- 
tative degree-degree correlations are introduced. The fi- 
nal fraction of nodes which ever hears a rumor (we call 
this the final size of a rumor), however, depends on the 
interplay between the model parameters and the under- 
lying topology. Our findings provide a better quantita- 
tive understanding of the complex interplay between the 
characteristics of a rumor, i.e. its spreading and cessa- 
tion rates, and the topology of the underlying interaction 
network that supports its propagation. They are rele- 
vant to a number of rumor-like processes such as chain 
emails, viral advertising and large-scale data dissemina- 
tion in computer and communication networks via gossip 
protocols. 

The rest of this paper is organized as follows. In Sec- 
tion II we describe our general rumor model. In section 
III a formulation of the model within the framework of 
Interactive Markov Chains is given, and the correspond- 
ing mean-field equations are derived. In section IV ana- 
lytical results are presented for the case of homogeneous 
networks, characterized by a constant degree and no de- 
gree correlations. This is followed by a numerical study 
of the model on ER random graphs and both uncorre- 
cted and assortatively correlated SF networks. We close 
this paper in section V with conclusions. 



II. A GENERAL RUMOR MODEL ON 
NETWORKS 

The model we shall consider is defined in the follow- 
ing way. A closed population consisting of N members is 
subdivided into three groups, those who have not heard 
the rumor, those who have heard it and wish to spread 
it, and those who have heard it but have ceased com- 



municating it. We call these groups ignorants, spreaders 
and stiflers, respectively. The rumor spreads by directed 
contact of the spreaders with others in the population, 
which takes place along the links of an undirected social 
interaction network G = (V,E), where V and E denote 
the vertices and the edges of the network, respectively. 
The outcome of contact of a specified spreader with (a) 
an ignorant is that the ignorant becomes an spreader at 
a rate A, and (b) another spreader or a stifler is that 
the initiating spreader becomes stifler at a rate a. The 
spreading process starts with one (or more) element (s) 
becoming informed of a rumor and terminates when no 
spreaders are left in the population. In the Daley-Kendall 
rumor model and its variants stifling is the only mech- 
anism that results in cessation of rumor spreading. In 
reality, however, cessation can occur also purely as a re- 
sult of spreaders forgetting to tell the rumor, or their 
disinclination to spread the rumor anymore. Following 
a suggestion in we take this important mechanism 
into account by assuming that individuals may also cease 
spreading a rumor spontaneously (i.e. without the need 
for any contact) at a rate 5. 



III. INTERACTIVE MARKOV CHAIN 
MEAN-FIELD EQUATIONS 

We can describe the dynamics of the above model on a 
network within the framework of the Interactive Markov 
Chains (IMC). The IMC was originally introduced in 
mathematical sociology, as a means for modeling social 
processes involving many interacting actors (or agents) 
[26j . More recently they have been applied to the dynam- 
ics of cascading failures in electrical power networks , 
and the spread of malicious software (malware) on the 
Internet [28|]. An IMC consists of ./V interacting nodes, 
with each node having a state that evolves in time ac- 
cording to an internal Markov chain. Unlike conventional 
Markov chains, however, the corresponding internal tran- 
sition probabilities depend not only on the current state 
of a node itself, but also on the states of all the nodes to 
which it is connected. The overall system evolves accord- 
ing to a global Markov Chain whose state space dimen- 
sion is the product of states describing each node. When 
dealing with large networks, the exponential growth in 
the state space renders direct numerical solution of the 
IMCs extremely difficult, and one has to resort to cither 
Monte Carlo simulations or approximate solutions. In 
the case of rumor model, each internal Markov chain can 
be in one of the three states: ignorant, spreader or stifler. 
For this case we derive below a set of coupled rate equa- 
tions which describe on a mean-field level the dynamics 
of the IMC. We note that a similar mean- field approach 
may also be applicable to other dynamical processes on 
networks which can be described within the IMC frame- 
work. 

Consider now a node j which is in the ignorant state at 
time t. We denote with p J u the probability that this node 
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stays in the ignorant state in the time interval [t, t + At], 
and with p 3 is = 1 — pi the probability that it makes a 
transition to the spreader state. It then follows that 



P i = (i - Atxy 



(i) 



where g = g(t) denotes the number of neighbors of node 
j which are in the spreader state at time t. In order to 
progress, we shall coarse-grain the micro-dynamics of the 
system by classifying nodes in our network according to 
their degree and taking statistical average of the above 
transition probability over degree classes. 

Assuming that a node j has k links, g can be considered 
as an stochastic variable which has the following binomial 
distribution: 



(2) 



where 9(k,t) is the probability at time t that an edge 
emanating from an ignorant node with k links points to 
a spreader node. This quantity can be written as 

k' k> 

(3) 

In this equation P(k'\k) is the degree-degree correla- 
tion function, P(sk'\ik) the conditional probability that 
a node with k' links is in the spreader state, given that 
it is connected to an ignorant node with degree k, and 
p s (k' , t) is the density of spreader nodes at time t which 
belong to connectivity class k. In the above equation 
the final approximation is obtained by ignoring dynamic 
correlations between the states of neighboring nodes. 

The transition probability pu(k,t) averaged over all 
possible values of g is then given by: 

p u (k,t) = ^( k )(i-\Atye(k,ty(i-e(k,t)) k -o 



3=0 



1- \AtJ2 p (k'\k)p s (k',t) 



(4) 



In a similar fashion we can derive an expression for 
the probability p ss {k,t) that a spreader node which 
has k links stays in this state in the interval [t,t + At}. 
In this case, however, we also need to compute the 
expected value of the number of stifler neighbors of the 
node at time t. Following steps similar to the previous 
paragraphs we obtain 



,(M) = (^-aAtY^P(k'\k){p s {k',t) + p r {k',t))^ 



x (1 - 6 At). 



(5) 



The corresponding probability for a transition from 
the spreader to the stifler state, p sr (k,t) is given by 

p sr (k,t) = 1 -p ss (k,t). 



The above transition probabilities can be used to set 
up a system of coupled stochastic equations for change in 
time of the population of ignorants, spreaders and stiflers 
within each connectivity class |30|. However, ignoring 
fluctuations around expected values we can also obtain a 
set of deterministic rate equations for the expected values 
of these quantities. In the limit At — > these equations 
are given by 



V(M) 



at 



'k\ P l {k,t)^p s {k',t)P{k l \k) (6) 



d P s (k,t) 
dt 



/cAp i (fc,t)^ / 9 s (fc',i)F(fc'|A:) 

k' 

kap s (k, t) t) + p r (k', t))P(k'\k) 



6p s (k,t) 



(7) 



d P r (k,t) 
dt 



= kap s {k, t)^2(p s (k', t) + p r {k', t))P(k'\k) + 
k> 

+ 5p s (k,t). (8) 



In the above equations p l (k, t), p s {k, t), and p r {k, t) are 
the densities at time t of, respectively, ignorants, spread- 
ers and stiflers in class k. These quantities satisfy the 
normalization condition p l (k, t) + p s (k, t) + p r (k, t) = 1. 

For future reference we note here that information on 
the underlying network is incorporated in the above equa- 
tions solely via the degree-degree correlation function. 
Thus in our analytical and numerical studies reported in 
the next section we do not need to generate any actual 
network. All that is required is either an analytical ex- 
pression for P(k'\k) or a numerical representation of this 
quantity. 



IV. RESULTS AND DISCUSSIONS 
A. Homogeneous networks 

In order to understand some basic features of our ru- 
mor model we first consider the case of homogeneous net- 
works, in which degree fluctuations are very small and 
there are no degree correlations. In this case the rumor 
equations become: 



dp*(t) 

dt 
dp s {t) 

dt 



d P r {t) 
dt 



= -\kp\t)p s (t) (9) 

= Xkp\t)p s {t) - akp s (t)(p s (t) + p r (t)) 

~ Sp s (t) (10) 

- akp s (t)(p s (t)+p r (t))+Sp s (t), (11) 



where k denotes the constant degree distribution of the 
network (or the average value for networks in which the 
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probability of finding a node with a different connectivity 
decays exponentially fast). 

The above system of equations can be integrated ana- 
lytically using an standard approach. In the infinite time 
limit, when s(oo) = 0, wc obtain the following transcen- 
dal equation for R = r(oo), the final fraction of nodes 
which ever hear the rumor (we call this the final rumor 
size) 



R = 1 



where 



(a + X)k 
S + ak 



(12) 



(13) 



Eq. (12) admits a non-zero solution only if e > 1. For 
S ^ this condition is fulfilled provided 

A: 



rfc > 1, 



(14) 



which is precisely the same threshold condition as found 
in the SIR model of epidemic spreading on homogeneous 
networks 0, |2t|. On the other hand, in the special 
case 5 = (i.e when the forgetting mechanism is absent) 
e = 1 + X/a > 1, and so Eq. (14) always admits a non- 
zero solution, in agreement with the result in p5|. 

The above result shows, however, that the presence of 
a forgetting mechanism results in the appearance of a fi- 
nite threshold in the rumor spreading rate below which 
rumors cannot spread in homogeneous networks. Fur- 
thermore, the value of the threshold is independent of a 
(i.e. the stifling mechanism), and is the same as that for 
the SIR model of epidemics on such networks. This result 
can be understood by noting that in the above equations 
the terms corresponding to the stifling mechanism are of 
second order in p s , while the terms corresponding to the 
forgetting mechanism are only of first order in this quan- 
tity. Thus in the initial state of spreading process, where 
p s w and p r xs 0, the effect of stifling is negligible rela- 
tive to that of forgetting, and the dynamics of the model 
reduces to that of the SIR model. 



B. Random graphs and uncorrelated scale- free 
networks 



Next we numerically investigate the dynamics of rumor 
spreading on complex networks. We consider first uncor- 
related networks, for which the degree-degree correlations 
can be written as 



P(k'\k) = q(k') 



k'P(k') 
~Jk~T' 



(15) 



where P(k') is the degree distribution and (k) is the av- 
erage degree. We consider here two classes of such net- 
works. The first class is the Erdos-Renyi random net- 
works, which for large N have a Poisson degree distribu- 
tion: 



P(k) 



(kl 
k\ 



(16) 



The above degree distribution peaks at an average value 
(k) and show small fluctuations around this value. The 
second class we consider are scale-free networks which are 
characterized by a highly right-skewed power law degree 
distribution: 



P{k) = 



Ak-~< 




otherwise. 



(17) 



In the above equation fc m i n is the minimum degree of 
the networks and A is a normalization constant. For 
2 < 7 < 3 the variance of the above degree distribution 
becomes infinite in the limit of infinite system size while 
the average degree distribution remains finite. We shall 
consider henceforth SF networks with 7 = 3. 

Our studies of uncorrelated networks were performed 
using the above forms of P(k) to represent ER and SF 
networks, respectively. The size of the networks consid- 
ered was TV = 10 5 N = 10 6 , and the average degree 
was fixed at (k) = 7. For each network considered we 
generated a sequence of N random integers distributed 
according to its degree distribution. The coupled set of 
differential equation (6-8) were then solved numerically 
using an standard finite difference scheme, and numerical 
convergence with respect to the step size was checked nu- 
merically. In the following and throughout the paper all 
calculations reported are performed by starting the ru- 
mor from a randomly chosen initial spreader, and averag- 
ing the results over 300 runs with different initial spread- 
ers. The calculations reported below were performed for 
networks consisting of N = 10 6 nodes. 

In our first set of calculations we set 5 = 1 and investi- 
gate the dynamics as a function of the rumor spreading 
rate A and the stifling rate a. First we focus on the im- 
pact of network topology on the final size of a rumor, R, 
which for inhomogeneous networks is obtained from 



fl = 5> r (M~), 



(18) 



where t^ denotes a sufficiently long time at which the 
spreading process has reached its steady state (i.e. no 
spreaders are left in the network). In Fig. 1 R corre- 
sponding to the ER network is plotted as a function of 
A, and for several different values of the stifling param- 
eter a. It can be seen that in this network R exhibits a 
critical threshold A c below which a rumor cannot spread 
in the network. Furthermore, just as in the case of ho- 
mogeneous networks, the value of the threshold does not 
depend on a, and is at A c = 0.12507. This value is in ex- 
cellent agreement with the analytical results for the SIR 
model on an infinite size ER network, which is given by 
A c = (k)/(k 2 ) = 0.125 We also verified numerically 
that the behavior of R in the vicinity of the critical point 
can be described with the form 



R~ A(X-X c f, 



(19) 



where (3=1, and A = A(a) is a smooth and monotoni- 
cally decreasing function of a. The results are shown in 
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FIG. 1: The final size of the rumor, R is shown as a function 
of the spreading rate A for the ER network of size 10 6 . The 
results are shown for several values of the stifling parameter 
a. 



FIG. 2: R is plotted as a function of A — A c for the ER network 
of size 10 6 , using different values of a. Solid lines show our 
numerical fits to the form R ~ (A — Ac)' 5 , with (3 = 1. 



Fig. 2 where R is plotted as function of A for a range of 
values of a, together with the corresponding fits. 

We further examined the above numerical findings by 
analytically solving Eqs. (6-8) in the vicinity of the crit- 
ical rumor threshold, to first order in a. Details of the 
calculations are given in the Appendix, and they confirm 
the above numerical results for A c , and the behavior of 
R in the vicinity of the critical rumor threshold. 

Next we turn to our results for the SF network. In Fig. 
3 results for R in this network are shown. In this case 
we also observe the presence of an a-independent rumor 
threshold, albeit for much smaller spreading rates than 
for the ER network. We have verified numerically that 
in this case the threshold is approached with zero slope, 
as can also be gleaned from Fig. 3. Since the value of the 
threshold is independent of a, we can use the well-known 
result for the SIR model (the a = case) to conclude 
that in the limit of infinite system size the threshold seen 
in the SF network will approach zero. It is therefore not 
an intrinsic property of rumor spreading on this network. 

In order to further analyze the behavior of R in SF 
networks, we have numerically fitted our results to the 
stretched exponential form, 

R ~ exp(— C/A), (20) 

with C depending only weakly on a. This form was found 
to describe the epidemic prevalence in both the SIS and 
the SIR model of epidemics 0, 0|. The results are 
displayed in Fig. 4, and they clearly indicate that the 
stretched exponential form also nicely describes the be- 
havior of R in our rumor model. This result provides 
further support for our conjecture that the general ru- 
mor model does not exhibit any threshold behavior on 
SF networks (at least in the limit of infinite systems size). 

In addition to investigating the impact of network 
topology on the steady-state properties of the model, it is 
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FIG. 3: The final size of the rumor, R is shown as a function 
of the spreading rate A for the SF network of size 10 6 . The 
results are shown for several values of the stifling parameter 
a. 



of great interest to understand how the time-dependent 
behavior of the model is affected by topology. In Figs. 5 
and 6 we display, as representative examples, time evo- 
lution of, respectively, the total fractions of stiflers and 
spreaders, in both networks for A = 1 and two sets of 
values of the cessation parameters: {8 = l,a = 0}, and 
{8 = Q,a = 1}. The first set of parameters corresponds 
to a spreading process in which cessation results purely 
from spontaneous forgetting of a rumor by spreaders, or 
their disinclination to spread the rumor any further. The 
second set corresponds to a scenario where individuals 
keep spreading the rumor until they become stiflers due 
to their contacts with other spreaders or stiflers in the 
network. As can be seen in Fig. 5, in the first scenario 
the initial spreading rate of a rumor on the SF network 
is much faster than on the ER network. In fact, we find 
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FIG. 4: R (in log scale) in the SF network of size 10 6 is 
plotted as a function of 1/A and several values of a. Solid 
lines are our numerical fits to the stretched exponential form 
R = B(a) exp(-C(a)/A). The network size is TV = 10 6 

that the time required for the rumor to reach 50% of 
nodes in the ER random graph is nearly twice as long as 
the corresponding time on the SF networks. This large 
difference in the spreading rate is due to the presence of 
highly connected nodes (hubs) in the SF network, whose 
presence greatly speeds up the spreading of a rumor. We 
note that in this scenario not only a rumor spreads ini- 
tially faster on SF networks, but it also reaches a higher 
fraction of nodes at the end of the spreading process. 

It can be seen from Figs, 5 and 6 that in the second 
spreading scenario (i.e. when stifling is the only mecha- 
nism for cessation) the initial spreading rate on the SF 
network is, once again, higher than on the ER network. 
However, unlike the previous situation, the ultimate size 
of the rumor is higher on the ER network. This behavior 
is due to the conflicting roles that hubs play when the 
stifling mechanism is switched on. Initially the presence 
of hubs speeds up the spreading but once they turn into 
stiflers they also effectively impede further spreading of 
the rumor. 



FIG. 5: Time evolution of the density of stiflers is shown 
on the ER (dashed lines) and the SF network (solid lines) 
when the dynamics starts with a single spreader node. Results 
are shown for two sets of values of the cessation parameters 
{a — 0,5 = 1} and {a — 1,8 — 0}. The network sizes are 
TV = 10 6 . 




time 



FIG. 6: Time evolution of the density of spreaders is shown for 
the same networks, model parameters and initial conditions 
as in Fig. 5 



C. Assortatively correlated scale- free networks 

Recent studies have revealed that social networks dis- 
play assortative degree correlations, implying that highly 
connected vertices preferably connect to vertices which 
are also highly connected [l(J- In order to study the im- 
pact of such correlations on the dynamics of our model, 
we make use of the following "local" ansatz for the 
degree-degree correlation function 

P(k'\k) = {l-f3)q(k')+(35 kk ,; (0 < (3 < 1). (21) 

This form has been used previously in recent studies of 
the SIR dynamics on correlated scale-free networks [3ll 



|32| , and allows us to study in a controlled way the impact 
of degree correlations on the spreading of rumor. 

Using the above degree-degree correlation function we 
numerically solved Eqs. (6-8) for a SF network charac- 
terized by 7 = 3 and < k >= 7. The network size was 
fixed at TV = 100, 000, and we used two values for the 
correlation parameter: /3 = 0.2 and (3 = 0.4. Fig. 7 dis- 
plays R as a function of A, and for a = 0.5,0.75, 1 (the 
value of 6 was fixed at 1). 

It can be seen that below A « 0.5 a rumor will reach a 
somewhat smaller fraction of nodes on the correlated net- 
works than on the uncorrelated ones. However for larger 
values of A this behavior reverses, and the final size of 
the rumor in assortatively correlated networks shows a 
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FIG. 7: The final size of the rumor is plotted as a function of 
A and for several values of a in the SF network of size 10 . 
Results are shown in the absence (solid lines) of assortativc 
degree-degree correlations and in the presence of such corre- 
lations. The correlation strengths used are (5 = 0.2 (short 
dashed lines) and (5 = 0.4 (long dashed lines). 



higher value than in the uncorrelated network. We thus 
conclude that the qualitative impact of degree correla- 
tions on the final size of a rumor depends very much 
on the rumor spreading rate. We also investigated the 
effect of assortative correlations on the dynamics (tem- 
poral behavior) of rumor spreading and found that such 
correlations slightly increase the initial rate of spreading. 



V. CONCLUSIONS 

In this paper we introduced a general model of rumor 
spreading on complex networks. Unlike previous rumor 
models, our model incorporates two distinct mechanisms 
that cause cessation of a rumor, stifling and forgetting. 
We used an Interactive Markov Chain formulation of the 
model to derive deterministic mean-field equations for 
the dynamics of the model on Markovian complex net- 
works. Using these equations we investigated analytically 
and numerically the behavior of the model on Erdos- 
Renyi random graphs and scale-free networks with ex- 
ponent 7 = 3. 

Our results show the presence of a critical threshold 
in the rumor spreading rate below which a rumor can- 
not spread in ER networks. The value of this threshold 
was found to be independent of the stifling mechanism, 
and to be the same as the critical infection rate of the 
SIR epidemic model. Such a threshold is also present in 
the finite-size SF networks we studied, albeit at a much 
smaller value. However in SF networks this threshold 
is reached with a zero slope and its value becomes van- 
ishingly small in the limit of infinite network size. We 
also found the initial rate of spreading of a rumor to be 
much higher on scale-free networks than on ER random 



graphs. An effect which is caused by the presence of hubs 
in these networks, which efficiently disseminate a rumor 
once they become informed. Our results show that SF 
networks are prone to the spreading of rumors, just as 
they are to the spreading of infections. 

Finally we used a local ansatz for the degree-degree 
correlation function in order to numerically investigate 
the impact of assortative degree correlations on the dy- 
namics of our model on SF networks. These correlations 
were found to increase slightly the initial rate of spread- 
ing in SF networks. However, their impact on the final 
fraction of nodes which hear a rumor depends very much 
on the rate of rumor spreading. 

The basic assumption underlying the mean-field equa- 
tions derived in this paper is that all vertices within a 
given degree class can be considered statistically equiva- 
lent. Therefore our results are not directly applicable to 
structured networks in which a distance or time ordering 
can be introduced, or there is a high level of clustering. 
We are currently working on more elaborate approxima- 
tions of our model which could take into account such 
properties, and in particular the presence of clustering, 
which is known to be an important feature of social net- 
works. 

Furthermore, in the present work we assumed the un- 
derlying network to be static, i.e. a time-independent 
network topology. In reality, however, many social and 
communication networks arc highly dynamic. An exam- 
ple of such time-dependent social networks is Internet 
chatrooms, where individuals continuously make new so- 
cial contacts and break old ones. Modeling the dynamics 
of rumor spreading on such dynamic networks is highly 
challenging, in particular when the time scale at which 
network topology changes becomes comparable with the 
time scale of the process dynamics. We also aim to tackle 
this highly interesting problem in future work. 

Acknowledgments 

M. N. acknowledges the Abdus Salam International 
Centre for Theoretical Physics (ICTP) for a visiting fel- 
lowship during which some amendments to this work 
were made. Y. M. is supported by MEC through the 
Ramon y Cajal Program. This work was supported 
by BT and the Spanish DGICYT Projects FIS2004- 
05073-C04-01. We thank Keith Briggs for reading the 
manuscript. 

APPENDIX A: CRITICAL RUMOR THRESHOLD 
IN ER RANDOM GRAPHS 

Eq. (6) can be integrated exactly to yield: 

p*(k,t)=p*(k,0)e- Xk ^\ (Al) 

where p l (k, 0) is the initial density of ignorant nodes 
with connectivity k, and we have introduced the aux- 
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iliary function 



#*) = 5>(*) / p^k^df^ / (( P s (k, t')))df 



In the above equation and hereafter we use the shorthand 
notation 



«o(*)»=5>(*)0(*) 



with 



(A3) 



(A4) 



In order to obtain an expression for the final size of the 
rumor, R, it is more convenient to work with <f>. As- 
suming an homogeneous initial distribution of ignorants, 
p l (k,0) — p l Q, we can obtain a differential equation for 
this quantity by multiplying Eq. (7) with q(k) and sum- 
ming over k. This yields after some elementary manipu- 
lations: 

d ± = l-« 6 -A*»)-ty 

1 - «c- Afc *(*')»l {(kp s (k,t')))dt',(A5) 



where, without loss of generality, we have also put p[, k 1. 

0, and Eq. (A5) 



In the limit t — > oo we have 



dt 



becomes: 

= l-^-***-))-*^ 

\_ {{e -Xk<Kt') )) \ ((kp°(k,t')))dt', 



(A6) 

where (p^ — lim^oo <p(t). 

For a — Eq. (A5) can be solved explicitly to obtain 
For a^Owe solve (A5) to leading order in a. 
Integrating Eq. (7) to zero order in a we obtain 



p s (k,t) = 1-e 



-\k<f> 



y* e «c*-f) 









-d I » : ' iJ -' ' I , - u "" > dt'+0{a). 

(A7) 

Close to the critical threshold both <p(t) and (p^ are 
small. Writing <p(t) — <poof{t), where f(t) is a finite 
function, and working to leading order in (poo, we obtain 



p s (k, t) ~ -SXk^ / e^-^fWdt' + O(^) + 0(a) 
Jo 

(A8) 



Inserting this in Eq. (|A6fl and expanding the exponential 
to the relevant order in 6^ we find 



= 0oo [A«fc» -S-X 2 ((k 2 ))(l/2 + a((k))I)cb 00 ] 

+ 0{a 2 )+0{ctt) (A9) 

where / is a finite and positive-defined integral. The non- 
trivial solution of this equation is given by: 



X((k))-S 



A 2 ((P))(i + aI«fc}»- 



(A10) 



Noting that ((k)) = (k 2 )/(k) and {(k 2 )) = (k 3 )/{k) we 
obtain: 



2(k)(^X-S) 
A 2 (fc3)(i + 2a /^ 



This yields a positive value for <poo provided that 



A > Jfe^ 
S ~ (fc 2 )' 



(All) 



(A12) 



Thus to leading order in a the critical rumor threshold is 
independent of this quantity, and is the same as for the 
SIR model. In particular, for 5=1 the critical rumor 
spreading threshold is given by A c = (k)/(k 2 ), and Eq. 
(A12) simplifies to: 



2(fc)(A-A c ) 
Vco A 2 (fc 3 )(A c + 2a/)' 

Finally, R is given by 

i? = 5>(fc)(l- e - Afe ^), 



(A13) 



(A14) 



and expanding the exponential in Eq. (A6) we obtain 

2(fc) 2 (A-A c ) 



R w ^P{k)Xk(p 



X(k 3 )(X c + 2aI)' 



(A15) 



which shows that R ~ (A — A c ) in the vicinity of the 
rumor threshold. 
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